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Spatio-temporal scaling for out-of-equilibrium relaxation dynamics of an elastic 
manifold in random media: crossover between diffusive and glassy regimes 
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We study relaxation dynamics of a three dimensional elastic manifold in random potential from 
a uniform initial condition by numerically solving the Langevin equation. We observe growth of 
roughness of the system up to larger wavelengths with time. We analyze structure factor in detail 
and find a compact scaling ansatz describing two distinct time regimes and crossover between them. 
We find short time regime corresponding to length scale smaller than the Larkin length L c is well 
described by the Larkin model which predicts a power law growth of domain size L(t). Longer time 
behavior exhibits the glassy regime with slower growth of L(t). 
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I. INTRODUCTION 

Fluctuations around macroscopically condensed states 
such as charge density waves (CDW) and flux line 
lattices in superconductors Q, often exhibit glassy dy- 
namics due to frustration between the elastic restoring 
forces originated from the stiffness of the ordered state 
and random pinning forces brought by impurities. Ther- 
mally activated process dominates the slow relaxation 
dynamics of such systems much as spinglasses (see, e. 
g., @,0,[E1) and super-cooled liquids [7[. 

An important basic problem in the studies of the glassy 
dynamics is the isothermal aging, i.e., relaxation at fixed 
temperature from initial states far from equilibrium 
Aging of elastic manifolds in random media has been 
studied theoretically by dynamical mean-field theories 
and numerical simulations [E EE [111, EE EE EI El- 
Presumably there exists a dynamical length scale L(t) 
which grows with time t such that the system is e quili - 
brated on the wavelengths smaller than L(t) [IH fl6l.Tl7| . 
In other word L(t) is the size of local equilibrium do- 
main. Roles of L(t) have been examined extensively in 
the context of aging of spin-glasses (see for instance, Ref. 
[3 HUES])- While L(t) grows algebraically without the 
random pinning forces, the frustration drastically slows 
it down. Typically one expects that the growth law be- 
comes logarithmic due to the energy barriers which grow 
with the length scale [H EE EE III EE EE H . The pur- 
pose of the present paper is to analyze aging of an elastic 
manifold in random media in terms of L(t) and investi- 
gate crossovers between the two characteristic regimes: 
so called the Larkin regime and glassy regime. 

A standard theoretical model to study the above men- 



tioned CDW- like systems is the elastic manifold model 
in random potential, e.g. the Fukuyama-Lee-Rice model 
[IE EH j given by the following Hamiltonian 



H 



dr 



± K \V0(r,t)\ 2 



h(r) cos (9(r,t) - /3(r)) 
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Physically the scalar field 9(r) at position r in the space 
are understood as the local fluctuation of the phase part 
of order parameter of the condensate, such as the CDW 
state. The first term with elastic constant k indicates 
the elastic deformation energy which is minimized when 
9(r, t) is spatially uniform in the absence of the second 
term, the random field energy. This sinusoidal random 
potential is a periodic function with respect to 9{r, t) 
reflecting the underlying periodicity of the condensate. 
Both amplitude h(r) and phase /3(r) of the random-field 
are quenched random variables with short-ranged spatial 
correlations. Hereafter (. . .) means a thermal average 
and 7TT means an average over the quenched randomness 
(samples). 

Let us recall here some basic static properties of the 
system which is understood better than the dynamical 
properties of our interests. It is believed that physical 
properties of this kind of systems, such as the roughness 
characterized by B(r) — {(9(r)—9(0)) 2 ), arc different 
on three distinct length scales [27|. First, in very short 
length regime, perturbative analysis of the effects of dis- 
order predict algebraic growth of the roughness with dis- 
tance r: -B(r) oc r 2 ^ with some roughness exponent 
Below four dimensions it is known that £ = (4 — d)/2 
[IE EE EH- Then the perturbative regime, which we 
call as the Larkin regime in the following, must be ter- 
minated at the so-called Larkin length L c over which the 
effect of randomness overcomes the elasticity. Then the 
so called random manifold regime begins [30l , l3ll . EE EH 
where many metastable states exist and the roughness of 
the system is characterized by a nontrivial roughness ex- 
ponent £rm- In much larger length scales, amplitude of 9 



2 



eventually grows beyond the period of the random poten- 
tial. If the periodicity is relevant, the system cannot gain 
more benefit of the potential energy at a cost of elastic en- 
ergy. Then the last regime called the Bragg glass regime 
begins. In three dimensions large wave length fluctuation 
is highly suppressed and B{r) is no longer expressed by 
algebraic functions r 2C - but by a certain logarithmic func- 
tion of the distance r [2?], HH, [35[ . Then it is said that 
the system is in the Bragg glass phase where the sys- 
tem has a quasi-long-range order (QLRO). In the case of 
single harmonic potential as the present model Eq. |T]), 
the end of the Larkin regime and the beginning of the 
Bragg glass regime coincide, i.e., the transient random 
manifold regime does not exist [13, [H] . This is because 
the potential has only single characteristic scale, that is a 
period 27r, and does not have another smaller scale which 
yields the upper bound of the Larkin regime such as short 
ranged correlation length of potential along ^-direction. 

In this paper, we study the out-of-equilibrium relax- 
ation dynamics of the elastic system in the periodic ran- 
dom potential, Eq. (fT]). This system shows different types 
of dynamics at different time scales. We show that each 
of these is related to equilibrium spatial property in the 
corresponding length scale. Although we discuss in par- 
ticular the case of a three dimensional system we may 
comment on systems in general d-dimensions. 

In the next section, we review the dynamics in the 
Larkin regime, which can be examined analytically. In 
the section III, numerical analysis of the structure factor 
is shown. In the section IV, we propose a scaling law 
which describes the crossover from the Larkin regime to 
the glassy regime. In the final section, we present con- 
clusions and remarks. 



II. POWER-LOW DOMAIN GROWTH IN 
LARKIN REGIME 

When phase fluctuation is very small, 8(r,t) -C 1, the 
Hamiltonian Eq. |T]) is reduced to the so-called Larkin 
model [28), which is exactly solvable. When the second 
term in the r. h. s. of Eq. |T]) is expanded by 8 up to 
the linear order, the disorder effect is represented by the 
quenched random force ?y(r) = — h(r) sin/3(r), which is 
supposed to be a random Gaussian number satisfying 



77(r) = 0, 77(r)77(r') = A h 6(r - r'), 



(2) 



with a finite variance A^. The overdamped Langevin 
equation is written as 



SH 
S9(r, t) 



*V 2 0(r,<)+r ? (r) + C(r,t), (3) 



where 7 is the friction coefficient and ((r, t) is a Gaussian 
white thermal noise satisfying 

(C(r, t)) = 0, (C(r, t)C(T',t')) = 2jT5(r - v')5(t - f), 

(4) 



T being the temperature of the heat bath. The average 
(• • • } is taken over independent noise realizations. 

The formal solution with the uniform initial condition, 
6(r, 0) = for all r, is expressed as 

0(q, *) = -/" dt'e-^-^h Wq ) + C(qj t >)\ , (5) 
7 Jo 

where the functions of q are the Fourier transformations 
of the corresponding functions of r. The structure factor, 
i.e., Fourier transform of the scalar correlation function 
is obtained as 



B(d,t) 



= T- 



dr(0(O,t)0(r,t))e^ = (\9(q,t)\ 2 ) 

1 2 
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1 - e~« L <*> 

2 



where 



L(t) = vWy. 



(0) 



(7) 



Equation ^ means that the two fluctuations owing to 
temperature and randomness are decoupled in the Larkin 
model. But the growth of the two are characterized by a 
single time-dependent length L(t) such that the system 
is equilibrated over wavelength shorter than L(t), i.e., 
B(q,t) ~ B(q,t -> 00) = T/nq 2 + A h /n 2 q 4 for qL(t) > 
1. 

From Eq. ([6|), we can evaluate the amplitude of phase 
fluctuation as 



a(t) ee (0(r,i) 2 >= / dqB(q,t) 
d T 



(8) 



Co 



d-2 



4 — d k 



where C\ and C2 are positive constants. The third term 
diverges in equilibrium (t — > 00, L(t) — > 00) below four 
dimensions. By writing a(t) oc L(t) 2C ' one can read 
off the the roughness exponent of the Larkin model as 
(Larkin = (4 — d)/2. In the long length (time) scale, 
L(t) 3> Lt ee ^/kT/A/j , the third term due to the 
quenched randomness is dominant and the second term 
due to the thermal fluctuation term can be ignored. 

Note that the growth law of L(t) given by Eq. (J7J) is 
the same as in the absence of the random potential, i. 
e., diffusive dynamics with dynamical exponent z = 2. 
It means that the quenched random potential does not 
bring pinning effects at the level of its linear approxima- 
tion. It is easy to see that this is the case at any higher 
levels of perturbative treatments of the random potential. 



III. GLASSY DYNAMICS IN NONLINEAR 
POTENTIAL 

The linear approximation adopted in the previous sec- 
tion breaks down for large <r(t), which occurs when L{t) 
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becomes as large as the so-called Larkin length, L c ; 



L c (A h ) 



C 2 



A h 



-l/(4-d) 



(10) 



This length is derived from the threshold condition; 
a(t) = (27rci) 2 , where cl is a constant, similar to the 
Lindemann's constant 37], indicating the border below 
which the nonlincarity of random potential can be ig- 
nored. The corresponding time scale is 



t c {A h ) = L c {A h f 1 /n oc A h 



2/(4-d) 



(11) 



Beyond the Larkin length nonlinearity of the random 
potential yields many metastable states. The potential 
barriers between them will significantly slow down the 
growth of <j(i). Hereafter we call this nonpurturbative 
regime just as the glassy regime which may corresponds 
to the random manifold or Bragg glass regimes. Now we 
study the growth of roughness in the glassy regime by 
numerical simulations. 



A. Simulations 

In practice we consider the lattice version of Eq. ([3]). 
The equation of motion for the phase at the lattice point 
Ti is 



-j Y, [»<(«) - W)\ 

jen.n. 

-/i iS in(0^(i)-A)+O(i), (12) 



where 



2YTS ij 6(t 



t>). 



(13) 



Hereafter we set the coupling constant J and the friction 
coefficient Y to unity. 

Here let us explain some details of the numerical sim- 
ulations. Phase variables 8i(t)'s are put on the cubic lat- 
tice with size N = L sys 3 = 128 3 and periodic boundary 
conditions are imposed in all directions. The phase /Vs 
are introduced as independent uniform random numbers 
between and 2ir. On the other hand the strength of the 
random field hi is set to a uniform value h for all sites 
[H, so that 



h 2 /2. 



(14) 



We investigate the relaxation dynamics at various val- 
ues of the random field h = 1.0 — 5.0 and temperatures 
T = 0.0,0.5,1.0 and 1.5 0. We numerically solve the 
Eq. (fT2"|) by the second order stochastic Runge-Kutta 
method [38||. In the initial state, 9i is set to for all i. 
Physical quantities are averaged over 128 runs at least, 
each of which has independent realizations of random 
phase Pi and thermal noise Q(t). 
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FIG. 1: Snapshots of the profile of the structure factor 
B(q,t). The time changes uniformly in logarithmic scale as 
t/0.34 = 4°, 4 1 , • ■ ■ , 4 6 from the bottom to the top. The sym- 
bols with error bars are data for L sys = 128 and dotted lines 
are for L sys = 256. The arrows connect the values in the 
Larkin model (bold curves) given by Eq. ((6]l and numerical 
data obtained in the sinusoidal model at the same time t. The 
vertical line shows q 2 = t^ 1 [47]. 



We define the lattice version of the structure factor 
Eq. j6]) as 



B(q,t)=N(\9(c l> t)\') 
with 0(q,t) = A^ 1 V . %(i)e iq - r * 0. 



(15) 



B. Results 

Figure Q] shows some examples of the profile of B(q, t). 
Components for all q are zero at t = and the struc- 
ture factor grow with time t. It can be seen that at 
larger q the amplitude saturates to ^-independent but q- 
dependent value while at smaller q the amplitude remains 
i-dependent but g-independent. The above observation 
suggests that there is indeed a dynamical length scale 
L(t) which grows with time t such that components sat- 
isfying qL(t) S> 1 become equilibrated: the system has 
become rough on short wavelengths but remains flat at 
larger wavelengths. 

A simplest scaling which connects the dynamical 
regime qL(t) <g; 1 and the static regime qL(t) ^> 1 may 

be 0, EJ IS EE H 



B(q,t) 



,-(d+2C) 



B(qL(t)) 



(16) 



where £ is the roughness exponent and the scaling func- 
tion behaves as B{x) ~ const, for x ^> 1 and B(x) ~ 
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x d+2( _ gy integrating over q one obtains the correspond- 
ing scaling form for 



<x(t) = 5>(q,i) = (0M) 2 > 

q/0 



as 



<r(i) oc L(t) 



2C 



(17) 



(18) 



(Note that the above scaling holds only if the system 
size L sys is sufficiently larger than L(t) for a given time 
t.) Indeed one can find easily that the Larkin model 
discussed in section |TT] satisfies these scalings exactly at 
T = 0. However, the real behavior will be more compli- 
cated even in the Larkin regime because roughness orig- 
inates not only from the quenched random field but also 
from the thermal noise at finite temperatures. Further- 
more, there will be a crossover from the Larkin regime 
to the glassy regime at L c . In section ITVl we perform a 
more elaborate scaling analysis taking into account these 
complications. 

The analytic solution of the Larkin model is also plot- 
ted in Fig.Q]for comparison. In very short time the struc- 
ture factors of the two models coincide. As time goes by, 
it becomes apparent that B(q,t) for the sinusoidal po- 
tential model Eq. (JT^J) docs not grow as fast as that of 
the Larkin model Eq. ©. As shown later, the time scale 
beyond which the equivalence breaks down is t c given by 
Eq. (jlip . Furthermore by a closer inspection it appears 
that the envelope function B(q,t — > oo), i.e., the equi- 
librium structure factor, changes from that of the Larkin 
model. The structure factor for small q parts seems to 
be slightly different from q~ A of the Larkin model. We 
regard these changes as the crossover from the Larkin 
regime to the glassy regime, which we analyze more care- 
fully in the section ITVl 

Figure [2] shows the time evolution of cr(t). In the 
Larkin regime L(t) <C L c or t <C t c , we expect cr(t) oc 
L(t) 2 < ~ t 2 ^ z (See Eq. ^TBJ)) with ( = 1/2 and z = 2. 
However, the data deviate from this behavior for t > t c . 
The crossover time increases as the strength of the ran- 
dom field h decreases. We consider this reflects crossover 
from the Larkin regime to the glassy regime. Indeed 
by simply scaling t by the anticipated crossover time 
t c oc h~ A given in Eq. (fTTj) . data of a(t) collapse onto 
a universal function for sufficiently large t. The growth 
of a(t) for t > t c is very slow presumably due to activated 
glassy dynamics. 



IV. SCALING OF STRUCTURE FACTOR 

As observed in the previous section, the crossover from 
the Larkin regime to the glassy regime can be described 
by a simple scaling in which the time t is scaled by the 
crossover time t c (Ah) given by Eq. pip corresponding to 
the Larkin length L c (A/j). Now we analyze the spatio- 
temporal scaling law of the structure factor B(q, t) itself, 
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FIG. 2: Time evolution of cr(t). The inset is a semilogarithmic 
plot of the same data. The solid line curves show the results 
of the Larkin model for h <C T. The horizontal dotted line 
indicates the equilibrium value a(t = oo) for h = 0. The lines 
noted with 2hcl and t c (h) indicate the crossover (see detail 
in [47]). 



which provides us more detailed information than the in- 
tegrated one a in Eq. (fT7|) . The basic idea is expressed 
by Eq. p^]) which connects the dynamic qL(t) <C 1 and 
static qL(t) 3> 1 regimes. However we need to take into 
account complications due to the roughness of different 
origins: i) thermal roughness with the roughness expo- 
nent (thermal = (2—d)/2, ii) roughness due to the random 
field in the Larkin regime with (Larkin = (4 — d)/2 and 
iii) roughness in the glassy regime which has different 
g-dependence with ii). 



A. Scaling ansatz 

We propose the following scaling ansatz. Most im- 
portantly the crossover from the Larkin to the glassy 
regime is taken into account by scaling the length (or the 
wave number) by the Larkin length L c given by Eq. (|10p 
and the time by the corresponding time scale t c given by 
Eq. pip . We propose that the structure factor takes the 
following form, 



B(q,t) = TL 



1 _ e -2Y(qL c )X(t/Ll) 



Y(qL c ) 

X _ e -Y{qL c )X(t/Ll) 

Y(qL c ) 



(19) 



This is an extended version of Eq. ©: the first and sec- 
ond terms describe the fluctuation due to thermal noise 
[45l | and quenched randomness, respectively. The dynam- 
ical and static regimes described in the simplest scaling 
Eq. (UHJ) correspond to XY <§C 1 and XY 3> 1 rcspec- 
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tively: in the dynamical regime XY <C 1 the struc- 
ture factor is ^-dependent but ^-independent while in the 
static regime XY 3> 1, it becomes i-independent but 
g-dependent. 

We suppose that the scaling functions X and Y take 
the following asymptotic forms in the Larkin regime, 



X(t) = 



t 



for J«l 



X g (t) for * > 1, 



Y(q) 



CSt. X <7 d / 2 +Cg 



for q 3> 1 
for 



(20) 
(21) 



The scaling function Y(q) with q = qL c describes the 
equilibrium structure factor: Y(q)~ 2 = q~ 4 in the Larkin 
regime q 3> 1 and <j~( d + 2 ?9) in the glassy regime q <C 
1. The exponent £ s is an unknown roughness exponent 
in the glassy regime which will be smaller than CLarkin- 
Particularly £ 9 will be zero if the system has quasi long 
range order. On the other hand, the scaling function 
X(t) with t = t/t c describes the growth law of L{t) in the 
Larkin regime i <^C 1 and the glassy regime t > 1. More 
precisely, the dynamical length L(t) can be estimated by 
solving 



Y(L{tr l L c )X{t/Ll) = 1. 



(22) 



The function X g (t) is an unknown increasing function 
of the scaled time t in the glassy regime which will be 
slower than any algebraic functions due to the anticipated 
activated dynamics [111, [l5| . 

For ic 1 and 5 » 1, the above scaling reproduces 
Eq. ^ in the Larkin regime. However it turns out that 
certain vertex corrections are needed in the coefficients 



T -> T(l 



a 7 



and 



A h (l 



,A h + ---) 



(23) 



(24) 



in analyzing the raw data. These coefficients appear 
when performing perturbation expansion of the random 
potential beyond the linear approximation in section [III 
In the following analysis we treated the first order correc- 
tion terms only and regarded them as fitting parameters. 



B. Numerical analysis 

Now let us examine the validity of the scaling law pre- 
sented above using our numerical data. Scaling functions 
X(i) and Y(q) are determined by least square fitting. 

First we perform fitting by fixing the temperature T. 
For example, Fig. [3] shows the result of the scaling for 
T = 0.5. We plotted F(q,t) as a function of XY where 



F(q,t) = - 



1 B(q,t) T 1-e 



-YX 
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FIG. 3: Result of the scaling at T = 0.5. F(q,t) is plotted as 
a function of X(t/t c )Y (qL c ) for h=1.0, 1.6, 2.0, 2.5, 2.7, 3.0, 
3.2, 3.5 and 4.0. The white curve indicates (1 - e~ XY )/XY. 
Here the coefficients are corrected with a T & = 0.027 as in 
Eq. ([23j) and a at = 0.22 as in Eq. (|24]). The rectangles in 
the inset show the data ranges of our simulations in t/t c and 
q 2 L 2 c for ft=1.0, 2.0, 3.0 and 4.0. The bold curve in the inset 



indicates X(t/t c )Y(qL c ) 
this curve is equilibrated. 



1. The space-time region above 



which leads (1 ■ 



-XY 



) / XY if the scaling law in Eq. |T 



is valid. The figure shows nice collapsing of data using 
proper scaling functions, X and Y. 

Next let us bring together data at different tempera- 
tures. To this end we note that the Larkin length can 
weakly depends on temperature, 



L c (A h ,T) = [A h /c L (T) 2 ] 



2l-l/(4-d) 



(25) 



LtA h LlA h 



Y 



The constant ct, (T) will be larger for higher temperatures 
because thermal fluctuations will mask the quenched ran- 
dom potential over short length scales. In practice we 
treated cl(T) as a fitting parameter. Then as shown 
in Fig. 21 scaling function Y for different temperatures 
can be laid on a universal curve independent of tempera- 
tures. The resultant correction factor cl(T) is shown in 
the inset of Fig. 2] 

The scaling function X(t) exhibits the crossover from 
the power-law domain growth to the glassy dynamics. 
Looking Fig. 0] carefully, the scaling function X is not 
universal, which is more apparent in lower temperature. 
Particularly the relaxation stops on the way at zero tem- 
perature. This indicates that the relaxation in the long 
time regime is thermally activated process. The system 
cannot escape from a metastable state without thermal 
assistance. 

The scaling function Y~ 2 represents the structure fac- 
tor in the equilibrium state (X — > 00). From the present 
scaling we obtain its shape even for the small wave num- 
bers q <C 1 where B(q,t) is still far from equilibrium. 
While the growth law X(t) turned out to depend on the 
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q 

FIG. 4: Scaling function X(t) and Y(q) for several fixed tem- 
peratures. The inset in the top graph indicates the correc- 
tion factor cl(T). Here we use cl(T = 0) by extrapolating 
from finite temperature data. More precisely speaking the 
horizontal axes in this plot, t and q, are scaled value using 
L c (A h ,T) = [hc L {l.0)/c L (T)]~ 2 instead of that in Eq. (l25i 

temperature T, we find that Y(t) is essentially indepen- 
dent of the temperature T. This means that spatial cor- 
relation function in equilibrium has universal form inde- 
pendent of both T and Ah- The long wave length behav- 
ior of Y seems to obey a power law; Y(q) ~ g-( d / 2 +C) _ 
The roughness exponent C is smaller than that in the 
Larkin regime and the value is consistent with the one 
for the random manifolds £ rm w 0.21 (46|. However it 
is more natural to expect that this agreement is a tran- 
sient behavior that £ is approaching to zero because the 
present single harmonic model is considered to take single 
crossover to the Bragg glass regime [13, |36| . 

V. SUMMARY AND DISCUSSIONS 

In this paper the relaxation dynamics of the three di- 
mensional clastic manifold in random potential has been 



studied. We especially focused on the crossover between 
the Larkin regime and the glassy regime, i.e., power- 
law domain growth and thermally activated relaxation. 
We proposed a new scaling method for the dynamical 
structure factor which encodes dynamical growth of the 
roughness of different origins and successfully applied it 
to explain the crossover between the Larkin and glassy 
regime. At a given temperature the structure factor 
out-of-equilibrium can be scaled by using the dynami- 
cal length L(t) and the Larkin length L c (Ah). Quite 
interestingly our analysis yields the structure factor in 
equilibrium, which is hard to observe in equilibrium sim- 
ulations. The temperature dependence can be also taken 
care only by introducing a correction factor c^(T) for the 
Larkin length. It turns out that the scaling function Y(q) 
are universal and independent of the temperatures. 

Although we analyzed the model with random po- 
tential that is periodic with respect to the deformation 
field, indication of QLRO was hardly observed. This 
will appear when the amplitude of phase fluctuation be- 
comes greater than the period of the random potential, 
cr{t) 3> (27r) 2 . The crossover region between the Larkin 
and Bragg glass regimes, however, persists for quite long 
time and has special importance in the dynamical as- 
pect. This is because glassy behavior becomes serious at 
the early stage of the crossover. From the obtained scal- 
ing function Y(q), we can roughly estimate the end of 
the Larkin regimes as cl — \/a(t c ) /2n « 0.1 in Eq. (fTO]) 
[47l ] . The growth rate of a(t) decreases quite quickly be- 
fore a(t) reaches (2-7r) 2 (see Fig. [2]). (In fact almost all of 
in our simulations are smaller than n and the system 
does not feel that the potential is periodic.) Therefore the 
early stage of the crossover, which has sufficiently long 
time range, hardly reflects the periodicity of the random 
potential and is similar to the crossover to the random 
manifold regime. 
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